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• The MSFC lunar impact monitoring program began in 
2006 in support of environment definition for the 
Constellation (return to Moon) program 

• Work continued by the Meteoroid Environment Office 
after Constellation cancellation 

• Lunar impact monitoring allows measurement of fluxes in a 
size range not easily observed 

• A paper published in Icarus reported on the first 5 
years of observations 

• Icarus: http://www.sciencedirect.com/science/article/pii/S0019103514002243 

• ArXiv: http://arxiv.org/abs/1404.6458 
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Observation Summary 




v> 



Morning observations 




330+ impacts 
since 2005 

Subset of 1 26 flashes on 
photometric nights to 2011 
141 hrs evening - 81 flashes 
126 hrs morning - 45 flashes 
Average: 2.1 hrs/flash 
evening/morning = 1.6:1 



Photometric error ~0.2 mag 
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Camera Field of View and Tracking 



Approximately 20 arcminutes 
horizontal, 3.8xl0 6 km 2 

Approximately lm effective focal length 
with Vi inch CCD 

Good compromise between collecting 
area and glare 

Use stars for photometric calibration 

Telescope mount with lunar rate (in RA 
and Dec) is helpful although manual 
corrections are needed 
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Automated Lunar and Meteor 
Observatory (MPC H58) 




• Telescopes 

14" (0.35m) 

Meade, Celestron 
Paramount (ME, MX) 

• Detectors 

Sony HAD EX - based video 
Gamma=0.45, man. gain, 
shutter off 




12/15/2006 
09:17:39.336 
33 ms 
m R = 7.4 
0.09 kg 

Geminid (35 km/s) 


11/17/2006 
10:56:34.820 
66 ms 

m R = 7.0 0 

0.03 kg 

Leonid (71 km/s) 


s 



11/03/2008 
00:11:06.144 
100 ms 
m R = 7.7 
0.1kg 

S. Taurid (27 km/s) 



04/22/2007 
03:12:24.372 
133 ms 
m R = 6.7 
0.08 kg 

Lyrid (49 km/s) 



T 
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Videos 
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Calibration: Magnitude Equation 


Parameters determined by observing stars with 
known magnitudes 

R = -2.5 log 10 (S) -k’X+T (B-V) + ZP 

R = Johnson-Cousins R magnitude 
k’ = extinction coefficient 
X = airmass (zenith = 1.0) 

T = color response correction term 

(B-V) = color index 

ZP = zero point for the night 

S = DN 1/0 45 if camera gamma set to 0.45 which 
extends dynamic range (faintest flash to saturation) 

DN = pixel value 0-255 
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Correction from HAD EX to R filter vs blackbody 

temperature 
R-EX replaces T(B-V) 



Theoretical peak flash temperature 2800K Nemtchinov et al. (1998) 
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Luminous energy from impact 

peak magnitude 


E lum=fx A ^ f X d 2 t Joules 
^ium = luminous energy 

AX = filter half power width, 1607 Angstroms for R 
/ = 2 for flashes near the lunar surface 
d = distance from Earth to the Moon 
t = exposure time, 0.01667 for a NTSC field 
f A = 10' 7 xlO + +zp R )/ 2.5 J cm -2 S -1 ^-1 

R = the R magnitude 

zp R = 0.555, photometric zero point for R (not the same 
as ZP in magnitude equation) from Bessell et al. (1998) 
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Mass of the impactor 

assuming impact speed (shower or sporadic 


Luminous efficiency 

r| = 1.5xl0~ 3 exp (-9.3 2 /v 2 ) 

v = impact speed in km/s 

Kinetic Energy 

KE = E lum I f] 

Mass 

M = 2 KE / v 2 
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Lgmlrt&uS Efficiency, 


Luminous Efficiency 




From Moser et al. (201 1 ) 
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Shower Correlation 
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Limiting Magnitude 
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Impact Kinetic Energy (joules) 


Impact Energies 



10 


10 


10 * \r 


10® k 


1 


: [ 



!! ! ! ! 


: ! : 


!! : : : 

!! ! I I 

II I II 

!! ! ! i 

* ■ ■ 

! i : 


i ih 


i 


\ + 


-i \± 


: : : : 


ts 


a: 

— * ■ 3 
yn ¥-s < 


& 

tu 


CD 

CJ 

1 £ 


TTT! !!! ! ! 


] i : : 


: i: : : 





i 


siin f 


120 ISO 240 

Solar Longitude (degrees) 


300 


360 


Red error bars - photometric uncertainty; Blue error bars - luminous efficiency uncertainty 
Squares indicate saturation 

The flux to a limiting energy of 1 .05 x 1 0 7 J is 1 .03 x 1 0' 7 km -2 hr 1 
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Limiting Mass 
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Meteoroid Masses 



Red error bars - photometric uncertainty; Blue error bars - range of reasonable luminous efficiencies 
Squares indicate saturation 


The flux to a limiting mass of 30 g is 6.14 x 10' 10 m' 2 y r 1 
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Comparison with Griin Flux 


• For our completion limit of 30g we saw 71 
impacts for a flux of 

6.14 x 10 10 m 2 yr 1 

• The Griin et al. (1985) flux above a mass of 30g is 

7.5 xlO 10 m^yr 1 
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Impact Flux at Earth Compared 
with Other Measurements 




fa- 

rt) 

0> 

> 

10 ? 

L. 

lit 

CL 

I0 e 

J= 

t 

ro 

10* 

UJ 

OJ 

10« 

JI 

■4—r 

J= 

■M- 

10 3 

5 

cn 

10 2 

_c 

■o 

10 1 

"5 

u- 

fa» 

10° 

J3 

E 

10- 1 

3 

c 

aj 8 

10' 2 

> 

m 

10 3 

3 

E 

3 

u 

1 

O 

i — i 

IQ-5 

1 


5atefl«te data {Brown at al. 2D02) 

8 BeVt fit line {Grown ct al. 2002) 

Small fireballs {Hallida? et aL 1996) 

Lunar oratemng {Werner et ai. 2002) 
LINEAR {Stewart 3001) 

LINEAR {Harms ZOQZj 
Spac-e^atch (Rabinowitz et al. 2000) 
neat (ftabmowiiz et al. 2QG0) 
Inlrasound/acousti-c bolides (Revel le 2001) 
Lunar Impacts, n -2e-l {Ortiz et al. 3006) 
Lunar Impacts, n*6e-3 {Orttz et al. 2006) 
Lunar Impacts, te^t (this work) 


Energy (kTTNT equivalent) 


After Brown et al. (2002) 
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Summary 


We have used a rigorous photometric procedure (observation of standards, 
color and extinction corrections, etc) to derive flash magnitudes 

- Brightest flashes are saturated; energy/mass under estimated 

Shower membership determined based on radiant visibility from impact 
location (zenith distance), time from maximum, and peak zhr (Figure of Merit 
in Suggs et al. 2014) - necessary for mass estimates 

Meteor showers are a significant contributor at cm sizes (>60%) - investigating 
radiant distribution as possible explanation for observed asymmetry 

Uncertainty in luminous efficiency dwarfs photometric errors - more work 
needed 

Flux results are consistent with other observational studies 
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Backup 
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Jack Schmitt/Apollo 17 
observation of lunar impact 



"NASA Apollo 17 transcript" discussion is given below (before descent to lunar surfac 

03 15 38 09 (mission elapsed time) 

(10 Dec 1972, 21:16:09 UT- possible Geminid) 

LMP Hey, I just saw a flash on the lunar surface! 

CC Oh, yes? 



Geminid visiblity 


LMP It was just out there north of Grimaldi [mare]. Just north of Grimaldi. You might see if you got anything on 
your seismometers, although a small impact probably would give a fair amount of visible light. 


CCOkay. We'll check. 


LMP It was a bright little flash right out there near that crater. See the [sharp rimed] crater right at the [north] edge 
of [the] Grimaldi [mare]? Then there is another one [i.e., sharp rimed crater] [directly] north of it [about 50km]- 
fairly sharp one north of it. [That] is where there was just a thin streak [pin prick] [flash?] of light. 


CC How about putting an X on the map where you saw it? 
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Meteor Shower and Sporadic Source Radiants 
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When We Observe 



• Initially, it was anytime the glare from the sunlit 
face did not completely wash out the earthshine 
face 

- Typically between 10% illuminated (crescent) and 50% 
(quarter) 

• Impact rate is higher during meteor showers and 
we are focusing on those now after 7 years of 
observing anytime 

• Observe from nautical twilight to moonset - 
evening 

• Observe from moonrise to nautical twilight - 
morning 

• Generate a schedule each year with dates, times, 
and shower visibilities 
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Equipment 



• Telescopes - 14 inch (0.35m), have also used 0.5m and 0.25m 

• Camera - B&W video l/2inch Sony HAD EX chip (Watec 902H2 
Ultimate is the most sensitive we have found) 

• Digitizer - preferably delivering Sony CODEC .AVI files if using 
LunarScan (Sony GV-D800, many Sony digital 8 camcorders, 
Canopus ADVC-110) 

— This gives 720x480 pixels x8 bits 

• Time encoder- GPS (Kiwi or lota) 

— Initially used WWV on audio channel with reduced accuracy 

• Windows PC with ~500Gb fast harddrive (to avoid dropped 
frames) 

— Firewire card for Sony or Canopus digitizers 
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Camera Settings 


• Manual gain control to do reliable photometry 

• Turn off automatic shutter control (ELC on 
Watec cameras) 

• No integration (Sense Up = off for StellaCam or 
MallinCams) 

• Best to use gamma = 0.45 to extend dynamic 
range at the expense of an extra calculation in 
the analysis (Gamma Lo for Watecs) 
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Optec 0.3x 
' focal reducer 


Watec 902H2 
. Ultimate 




Operator position 


South Dome Video 


Tower Video 
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Photometric Calibration 


• Use "all sky" photometry 

- Require standard stars with various colors at 
various airmasses 

• Calibration using earthshine is a bad idea 

- Brightness changes with terrestrial weather 

- Color changes with terrestrial weather 

- Extended source vs point source difficulties 

• Color correction between filtered magnitude 
of standards and color of flash is important 
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Comparison Stars 



• Stars will pass through the field of view during 
observations, but 

- you don't typically know the R magnitude 

- they are seldom in the FOV at the time of the flash 

- this means you must do "all sky" photometry rather than 
"differential" (i.e. must account for extinction as a function of 
airmass as well as zero point) 

- flat field must be very good because vignetting is worse near the edge of the 
FOV where the field stars will be seen, especially with focal reduction 

• Observe some "standards" at various airmasses (1 and 2 - 
3) after evening observations and before morning ones 

• Build a standards list using SIMBAD for stars that are 
bright enough but don't saturate the system (8 - 9 R mag 

for 14in) that pass through the zenith simbad.u-strasbg.fr/simbad/ 

- Must have published R and B-V mag and not be a variable or multiple 

R = -2.5 log® -k'X+T[B-V)+ ZP 
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Comparison Stars 


Filters and Photometric Calibration^^^ 


• Use the camera unfiltered to give maximum 
sensitivity 

- Wider spectral response 

• near infrared where the flash is brightest 

• blue and green where earthshine is brightest (to see features) 

• Calibration should be done with R magnitudes of 
comparison stars 

— Peak sensitivity of HAD EX and R filter is at the same 
wavelength but width is very different 

- Need the color term T (B-V) = EX-R in the magnitude 
equation 

R = -2.5 logfSj -k'X+T {B-V) + ZP 
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Software we have used 



• WinDV for recording windv.mourek.cz 

• LunarScan detection software (Gural will discuss) 

www.lunarimpacts.com/lunarimpacts.htm 

• VirtualDub for slicing out relevant sections of video and 
converting to "Old AVI" for reading into Limovie 

www.virtualdub.org/download.html 

• Limovie for checking photometry of flashes and 

calibration stars www005.upp.so-net.ne.jp/k_miyash/occ02/limovie_en.html 

• MaximDL can convert video segments to FITS 

- Don't use the aperture photometry tool until after each pixel 
is gamma corrected by 5 = DN 1/045 if camera gamma set to 
0.45 

• Python and Pyraf may be used for aperture photometry 

www.stsci.edu/institute/software_hardware/pyraf/current/download 
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Important Points 



• Flux determination requires a measurement of 
the number of hours of observation to a 
particular limiting magnitude 

- Do not count cloudy hours 

- Cumulative peak magnitude diagram is useful in 
determining limiting magnitude 

• We use peak magnitude rather than a time 
integrated magnitude 

- Later phases of an impact "light curve" are dominated 
by cooling of the ejecta and crater - relation to impact 
energy is contaminated by regolith properties 

- How long the flash is visible depends on variables 
such as atmospheric transparency and earthshine 
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